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SUMMARY 


A new, fast, non-iterative version of the “Wall Pressure Signature 
Method" is described and used to determine blockage and angle-of-at tack wind 
tunnel corrections for h igh 1 y-powe red jet-flap models. The correction method 
is complemented by the application of tangential blowing at the tunnel floor 
to suppress flow breakdown there, using feedback from measured floor pressures 
This tangential blowing technique was substantiated by subsequent flow invest- 
igations using an LV. 


The basic tests on an unswept, knee-blown, jet flapped wing were supple- 
mented to include the effects of s la t- remova 1 , sweep and the addition of 
unflapped tips. Cy values were varied from 0 to 10 and free-air 1 s in excess 
of 18 were measured in some cases. Application of the new methods yielded 
corrected data which agreed with correspond i ng large tunnel “free air” results 
to within the limits of experimental accuracy in almost all cases. A program 
listing is provided, with sample cases. 


The present report is the first of two parts: Part Two describes an 

extension to include jet- in-crossf low effects. A copy of the present report 
is retained in the Lockheed-Georg i a Company Engineering Report Files. The 
identifying number is !.G£»ER0166. 
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1.0 INTRODUCTION 


1 . 1 Background 

In any wind tunnel test, the basic requirement is to create a flow field 
around a test model which properly represents either free air conditions or, 
on occasion, the condition of flight near the ground. For conventional models, 
nominal tunnel velocity must be corrected in magnitude and direction to com- 
pensate for the presence of the tunnel walls. For V/STOL models these 
corrections are likely to be large enough to require special correction 
methods and the further complication arises that separation may be induced on 
a tunnel surface. If an in-ground condition is to be simulated the relative 
ground motion must also be considered: in flight, this motion will usually 
reduce the extent of a ground separation (if present) but will not necess- 
arily eliminate i . 

Within the above terms-of- reference , three distinct but related test 
needs may be identified: 

(a) the need for improved correction methods, particularly for 
blockage effects, including the effects of highly three 
dimensional powered flows. 

(b) the need to understand and either correct for or remove the 
effects of tunnel flow breakdown during tests to determine free 
air da t a . 

and (c) The need firstly to understand and then to properly simulate 
the effects of ground motion during ground effects testing. 

References 1 through 10 represent some ten year's work at Lockheed- 
Georgia on the above questions. As a result of this and the present work, 
the flow physics is now well understood and practical solutions are almost 
complete. To place the present work in perspective a review is presented 
below covering blockage experiments, software development, ang 1 e-of-attack 
correction and ground or tunnel floor separation phenomena as studied at 
Lockheed-Georg i a during the 1970's. 

Experiments on Wind Tunnel Blockage 

The history of wa 1 1 -pressure based tunnel blockage correction research 
at Lockheed is represented chrono 1 og i ca 1 1 y by References 5 through 10, or 
parts of these. 

When conducting an investigation of ground effects on a knee-blown 
flap model (Ref 2) a substantial static pressure drop was noticed between 
the test section entry and the tunnel breather slots at the test section 
exit. The calibrated velocity, at the test section entry, was evidently 
significantly below the effective value at the model implying that the 
conventionally calculated model coefficients were too high. An obvious 'fix' 
was to define a model station reference static pressure equal to the mean 
of the test section entry and exit values. 


This approach was applied to pressure data from rsew tests on a knee- 
blown flap model in the 30" x k2" tunnel and comparisons were made with 
datum tests in the Lockheed 16£' x 23i* tunnel (Ref 5). In the absence of 
balance data, CL-valjes were estimated from pressure integrations. Only a 
basic, straight wingeo, slats-on conf i gu ra t i on was tested. These pilot 
studies showed significantly improved C[_ cor re 1 a t i ons , between tunneis, 
when the new reference static correction procedure was employed. 

A fuselage containing a three-component balance and optional, unflapped 
wing tip extensions were added to the knee-blown flap model for the next test 
series (Ref 6). Test conducted in the 30*' x ^4 2 1 * tunnel and datum tests in 
the NASA/AAMRDL 7 ' x 10 1 tunnel included wake flow as well as balance measure- 
ments. With the slats fitted, the flow measurements showed little wake dis- 
tortion, relative to a corrected mainstream vector, and good force and moment 
correlations were obtained. However, with the slats removed the drag behavior 
in the small tunnel was totally different from that in the large tunnel, 
though the lift performance was comparable. Slats-off flow data were noL 
taken but analysis of the drag data suggested that flow breakdown in the 
short test section of the small tunnel interacted in some way with the 
separated main wing flow and caused the jet sheet to separate prematurely 
from the flap upper surface. In addition to this problem, it was recognized 
in Reference 6 that the revised reference static method responded primarily 
to wake blockage and was inherently incapable of responding to solid or 
separation -bubble-induced blockage. 

Reference 7 describes early Locki eed- f ended work on what has become 
known as the “wall pressure signature method." As the name indicates, a 
series of pressures along the test section length is used to characterize 
the tunnel flow. Analysis of this "signature" yields not only individual 
estimates of solid/bubble and wake blockage but also corresponding axial 
velocity interference increments anywhere in the test section. The feas- 
ibility of the approach was established by means of tests on normal flat 
plates of various sizes tested in the Lockheed l6i' x 23i ( tunnel. The 
data of Reference 7 were analyzed entirely by ‘hand 1 methods, using look-up 
charts: it was a considerable time before the cor res pond i ng computerized 

version was ready for ‘production’ use. 

From the work of References 6 and 7 it became clear that the L-foot 
test section length of the 30" x k2" tunnel was insufficient. "he tunne' 
test section was therefore reworked to 7~foot total length. Rows of 
permanent wall pressure orifices wsie edded. 

Reference 8 closely parallels Reference 6 but describes tests on a 
swept wing variant of the knee-b lewn- f lao model. The straight winged model 
was retested, in the longer test section, and the 'drag flip back' anomaly 
disappeared. The correlations for the straight wing improved and those f or 
the swept wing were good for at tached- f low cases. Wall pressure signatures 
were measured but were not used for correction purposes. Nonetheless, they 
gave important insight into tunnel interference and tunnel flow breakdown 
phenomena . 
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Other, Lockheed I RAD-sponsored , tests at this time included work on 
spheres of two sizes in two tunnels and on flat plate wings of four sizes 
tested in the (now) 30" x *0" tunnel- Automation of the wall pressure sign- 
ature method was completed in 1977 and its usefulness in application to auto- 
mobile testing in the 1 64 1 x 23i’ tunnel was becoming appreciated. However, 
it could be used only off-line because its operation was somewhat slow. 

Reference 9 collects together most of the previous data and analyzes 
it using the automated program, which it also documents. Data for normal 
flat plates, spheres, and idealized automobile, flat-plate wings and the 
unswept knee-blown flap model are all included. 

Software Development 

The initial objective of the computer program is to locate a source- 
sink pair, representing solid/bubble blockage and a wake-sou rce , all on the 
tunnel axis, and determine their strengths so as to provide the best curve 
fit to the observed wall pressure signatures. This is essentially an in- 
verse problem and the solution must be found iteratively with regard to the 
source and sink locations. A developed version of the previous look-up charts 
(Reference 7) is used, in tabular form, during this iteration. Having 
solved this inverse problem, the de t€ rm i nat i on of tunnel interference effects 
is s tra i ghtforward . 

The period from 197 d to 1978 saw substantial improvements in program 
capability with regard to increased robustness and reduced run time. It 
was found that a good deal of data reviewing is required to reject 'bad' 
points, to interpret unusually shaped signatures properly and to achieve 
the best theoretical match to observed data. The earliest program ran 
about 30-seconds per data point, which is totally unacceptable for on-line 
use. The Reference 9 program requires about 3~seconds on a minicomputer and 
is much more robust than the early programs. A practical limit appears to 
have been reached in development of the method in its iterative form. 

Reference 10 describes the most recent Lockheed research on the wall 
pressure signature method. An alternative approach to the iterative method 
is introduced in which multiple sources or sinks are employed at fixed 
positions. This method avoids iteration and a constant influence matrix 
may be used. A 1 eas t -squares fit to the wall pressure signature may be 
achieved, when using the new program, by choosing fewer singularities than 
pressure data points. The direct method is an order-of-magn i tude faster 
than the best iterative program. It can also accommodate unusually-shaped 
wall signatures for which the previous method must make approx imat ions . 

Angle -of- Attack Corrections 

The sensitivity to ang 1 e-of-at tack correction is either zero or weak 
in most of the correlations described above. It has been found sufficient 
to employ the methods of Williams and Butler (Reference ll) for the Dowered 
model tests or the classical, Glauert correction as quoted in Reference 12 
in other cases. However as is pointed out in Reference 10 t the development 
of a wall pressure signature method for ang 1 e-of-at tack is desirable to 
afford consistency with the blockage corrections. 

Referenc "9 describes intitial studies of ang 1 e-of-at tack correction 
by the wall pres^ire signature method. The general feasibility is estab- 
lished and a number of sensitivity studies are described. However, only 
limited examples are quoted which involve test data. 
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Ground or Tunnel-Floor Separation 

References 1 through 6 deal predominantly with ground simulation in 
the wind tunnel. It is clear the the most realistic simulation should 
include the ground motion, using a moving belt or some alternative means 
of controlling the flow immediately above the tunnel floor. It is shown 
in Reference 3 that tangential blowing along the qround, from just ahead 
of the model, may be used successfully to simulate a moving belt. The 
criterion for blowing quantity, based upon the physics of the flow, is 
that skin friction at the ground shall be positive or zero everywhere. 
Reference 6 describes the development of a ground blowing system which 
employs feedback from ground skin friction sensors to determine the level 
of tangential blowing. 

A blowing rig designed for ground simulation by tangential blowing 
(e.g. Ref. U) may also be used to control tunnel flow breakdown. However, 
there is an important distinction between the two applications. It is 
shown in Reference 3 that, even with a movinq qround, a spanwise vortex 
may be trapped between the wing and the ground. The appearance of a floor 
vortex during center-tunnel testing heralds the onset of tunnel flow break 
down and can never be a "correct" flow condition. We shall see later in 
this report that such a vortex can distort the flow seriously in the 
vicinity of the model and render the data uncor rec tab I e . When used during 
center tunnel testing, we shall see that ground blowing should be used to 
destroy the floor vortex, if it occurs. 


1.2 Some Theoretical Considerations 

Seleo^zcn o r Flow Model 

It is possible, if only in principle, to exploit the non-iterative, 
matrix approach described above by defining three-dimensional arrays of 
sources and vortices clustered in the vicinity of the model and its wake 
and solving for boundary conditions derived from wal I pressures. (In the 
present work, the normal velocity condition is satisfied by using an approp- 
riate image system.) However, such an approach would almost certainly 
encounter matrix conditioning problems. 

The number of unknown source or vortex strengths is reduced greatly if 
knowledge of the model geometry - location, wing sweep, angle-of-attack 
etc - is exploited. This relieves the matrix conditioning problem signifi- 
cantly though, as indicated in Reference 10. some difficulties remain. The 
problem becomes more one of limiting the number of influence matrices which 
must be held ready for use. 

Even after reducing the number of s I ngu l a r i t i es , there are constraints 
on their geometry which must be recognized. For example, if measured axial 
velocity at the tunnel wall mid height is used as a boundary condition, the 
strength of a vortex at mid height can not be determined because it cannot 
affect the boundary points' axial velocity. Sources Qj and Q 2 or vortices 
f * . p, placed at altitudes ih cannot be resolved separately because 

the boundary velocity depends only upon (Qj + Q 2 ) and (r 1 - r 2 ) respectively 
For this reason the inclusion of a vortex on the centerline or the inclusion 
of sources or vortices equally spared above and below the center I . ne results 
in a singular influence matrix for mid-wall orifice locations. These 
considerations suggest some necessary rules for valid singularity arrangements. 


t.* J l q ly w 




Other rules are probably needed to complete a set which is also sufficient: 
further work is needed to identify these. 

Uniqueness of the Interference Flow Field 

The constraints above reduce the permissible number of singularities, 
they restrict their location to the general area of the model and it’s wake, 
they introduce some geometric properties related to the model itself and they 
introduce certain restrictions intended to avoid singular influence matrices. 
Even within these constraints a considerable number of possible arrangements 
cf s I ngular i t i es remains, particularly with regard to their number and 
spacing. The details of the conf i gurat ion selected will affect the sing- 
ularity strengths but the implications in relation to the calculated inter- 
ference velocities are not immediately clear. 

Experience suggests that the interference flow field may be relatively 
insensitive to the fine details of the flow model. For example, a study 
is reported in Reference 10 in which the original source-source-sink, 
variable geometry formulation of the present problem was set up in non-linear 
equation form. A range of solutions was found, with widely varying geometry, 
and an interference velocity profile along the tunnel axis was calculated for 
each. Though the interference curve certainly was not unique, the spread 
between the individual solutions was acceptable in engineering terms. 

Interpretation of the Inter fere>ice Flow Field 

Having solved the inverse problem, as indicated previously, and having 
defined interference velocities at locations of interest on the model, what 
remains is to determine their effects. This subject is discussed in detail 
in Reference 10. if the maximum benefit is derived from the wall pressure 
signature method, wind tunnel models may be sufficiently large that inter- 
ference qradient corrections need to be cons i de r ed . If the pressure 
gradients are nearly constant it usually will be possible to use standard 
gradient correction ("ouoyancy") methods. If only surface pressure measure- 
ments are to be corrected, a "local mainstream" concept has been found to 
be effective in correcting for blockage (Ref 10). Beyond these, a method 
must be found for distributing the forces over the model so that moment 
corrections, in particular, may be made on the basis of the local conditions 
which apply to individual model components. Though an experimental approach 
is a candidate for this, and is used occas i ona 1 1 y , a better choice is 
probably a simple analytical model of the conf i gu ra t i on concerned. 

Once a h i gh- i nduced-g rad i ent field has been defined - by whatever 
method - it is highly desirable to seek out and exploit such flow models 
as are available for the conf i gurat ion concerned. A close interface with 
the "customer” is lively to be very beneficial in this regard. 

Impingement Cases for Powered Flows 

Even if the vortex which occurs ahead of floor impingement is removed by 
floor blowing boundary layer control, as described in sub-section 1.1, there 
may Still b t. sufficient flow distortion to make data correction difficult, 
if not impossible. However, with the floor vortex removed, there is at least 
a reasonable change of defining the interference flow field over the model 
vo l ume . 
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The calculation of the interference flow field for an impinging jet 
include^ the determination of the effects of truncation as well as the effects 
of images The vortex pair which might represent an impinging jet-in- 
crossflow bends sharply at the floor. To complete the interference calculation 
a contribution from the "missing" part of the plume downstream of this must 
be added to the image effects corresponding to the section of plume within 
the tunnel For flow continuity a further source effect, at the tunnel flo ^ r 
itself, may be needed to provide an appropriate envelope around the impinged 
jet fluid there. 

With the interference flow field defined, a final consideration concerns 
jet path distortion. To first order, this will be a jet velocity ratio 
effect wh.ch should be adequately accommodated when the corrected mainstream 
velocity is defined at the model location. As the plume of an impinging jet 
is likely to be aerodynamicalty "stiff" the distortion due to gradients bet- 
ween the model and the tunnel floor are likely to be insignificant to within 
a short distance from impingement. 


1.3 Layout of the Present Report 

This is the first volume of a two-part report. The present volume deals 
with conventional, winged configurations and includes computer program list- 
ings relevant to the baseline, wall pressure signature program. It should 
be noted that the baseline program is not restricted to unpowered cases: it 

will accommodate jet-flapped configurations, for example. Volume 11 dea s 
with the special topic of jet-in-crossflow modeling, as it affects wall 
pressure signature analyses. 

Section 2 of the present report comprises a description of the new, 
direct version of the wall pressure signature method. This repeats some 
Reference 10 material, but this is included for ready reference in connection 
with the corresponding program listings. 

Test hardware for recent knee-b 1 own- f 1 ap (KBF) model tests is described 
in Section 3' j e t - i n-c ross f I ow hardware details may be found in Volume II. 

The application of tangential blowing at the tunnel floor in KBF tests is 
described in Section h. 

Most of Section 5 comprises a presentation of results for several 
configurations of the knee-blown flap model and shows the correlations 
between 30" x ky' tunnel corrected data and cons t ra i nt - f ree data. The mam 
text of this report is completed by Discussion, Conclusion and References 
in Sections 6, 7 and 8, respectively. The Appendices include the appropriate 
p rog ram listings, user guides and data tables. 

The present report is intended to complement and update Reference 9 
which includes more detail on how the basic wall pressure signature method 
works together with practical details concerning its implementation. 


2.0 THE GENERALIZED METHOD FOR WALL PRESSURE SIGNATURE ANALYSIS 


2 . 1 Rev i ew 


Reference 7 includes the original formulation of the problem of 
determining wind tunnel blockage via the solution of an inverse problem, 
starting with measured wall pressures. The general approach is to find the 
strengths of an array of line sources and sinks, located on the horizontal 
center plane of the tunnel, which when acting with the appropriate wall image 
set, produces the observed wall pressures. Having solved this inverse pro- 
blern, tunnel blockage is determined by considering the image set acting alone. 
This approach is retained for the present work. An iterative solution which 
has been the standard approach to date is described, in its most developed 
form, in Reference 10. The more recent generalized, or matrix solution is 
also described in Reference 10 and sensitivity studies, to source or vortex 
span, phase etc. are also described. 

The algorithms for influence coefficient calculations were relatively 
s t ra ightforward in the Reference 7 and 10 programs, since only spanwise line 
sources were involved. However the geometric requirements for swept wing 
and for jet- i n-crossf low models are more demanding and a generalized, skewed, 

I i ne-s i ngu 1 ar i ty algorithm has been prepared. The formulat'on for sources, 
horseshoe vortices and doublets and the corresponding algorithms are docu- 
mented in Reference lA. 

In the sub-section which follows, some of the more important character- 
istics of the matrix approach will be reviewed. Some recent findings con- 
cerning the choice of pressure sensing points will be discussed in sub- 
section 2.3* The effects of model offset and sweep, on measured signatures 
are reviewed in subsection 2 . *4 and a 1 eas t -squa res formulation of the basic 
problem is given in sub-section 2.5. The section concludes with a mathe- 
matical description of the generalized method. 

2.2 Properties of the Influence Matrices and Their Inverses. 

Figure 2.1 shows influence matrices for five-element line-source and 
five-element horseshoe vortex systems. The source matrices are, in fact, 
the sum of two others, corresponding to the direct influence of the line 
sources (an antisymmetric matrix) and the influence of matching, but opposite- 
sign, sources situated far downstream which are needed to satisfy continuity. 
Every element in the downstream source matrix equals 0.5* Each of the 
constituent matrices is singular, but their sum is not. Inspection of the 
tunnel Moor and roof source coefficients (Figure 2.1) shows that, to avoid 
repeated rows in the influence matrix (which would make it singular) mean 
values of supervelocity increment must be determined from floor and roof 
orifices having the same x-location. However, sidewall data will generally 
be used for blockage estimation. 

The vortex influence coefficients include vertical velocity components 
at the tunnel sidewalls, as denoted by arrows in the upper right portion of 
Figure 2.1. Though these components could, in principle, be measured and 
used to determine vortex strengths this is less practical than measuring 
solely static pressures. However, we shall see later that these velocities 
may influence pressures s i gn i f i cant 1 y and hence may represent a I i f t -dcpcnden t 
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interference upon the blockage signature in some cases. The roof end 
floor vortex coefficients are of opposite sign, at a given x-location. 

When solving for lift interference, differences must be taken between 
supervelocity data determined at corresponding roof and floor orifice 
loca t i ons . 

Figure 2.2 shows a wail influence matrix for sources (upper left] I and a 
roof/r'loor influence matrix for vortices (upper right) together with their 
respective inverses, below them. In both cases the inverses mciuae anei- 
natinq-sign elements, indicating that the influence matrices are i 1 
conditioned. Though it has been demonstrated that correct singularity 
strengths are returned from computer-generated wall pressure signatures, 
may be anticipated that, for -noisy 1 , real data, oscillating S ’Tda^ con- 
strengths will be returned. Application of the method to tunnel data con 

firms this (see Ref. 10). 

To complete a tunnel interference calculation, the source or vortex 
effects at the tunnel centerline are determined, with the centra system 
removed This step may be combined with the previous one by multiplying 
the°cen ter -tunnel interference matrix by the inverse matrix already deter- 
mined. The product matrices are shown in the lower part of Figure 2.2 As 
before the elements have alternating signs. Nonetheless it is found t 
smooth’ interference distributions are generally obtained from experimental 

data . 

Figure 2.3 shows results from pilot tests on an interim program, 
designated -'MATCH", which employs the new matrix method. Correspon i ng 
results are also shown using the previous iterative program. The wal 
pressure signature fitted by "MATCH" passes through every experimental 
point: the iteratively obtained signature must approximate because « _ 

fewe- degrees of freedom. Though the source-sink geometries differ co 
siderably, the two methods predict remarkably similar distributions of 
interference velocity. 

2.3 Geometrical Considerations 
Si yiaul ivitu Svacinq and Location 

In ea‘rly P studies. solutions were obtained using arrays like those 
shown in Figure 2.1. Though good interference prediction £ as . p ^ s ® 

(Fig 2.3) wildly oscillating singularity strengths were obtained which 
were obviously unrelated to the flow physics. Closing up the arra * s ^ 
placing them around the model location would, in principle, rel.e 
problem but in practice did not because the matrix became increasingly ill- 
conditioned. It is evident from Figure 2.3 that a satisfactory solut on 
is obtained with a reduced number of singularities, provided that their 
placement recognizes the model and flow geometry appropriately. To satis y 
the greater number of boundary condi tion> a least-squares approach is the 

fore required. 

In addition to the downstream sink, matching each source in the test 
n r : x * sinqle upstream source is also provided, explicitly* to 

h e. P , » « hi.« . 

.. . 


allow the overal I signature 10 Mini ' 

better match to the experimental upstream asymptote. 
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Singularity Span 

it has been found that the present, generalized method is fairly 
forgiving with regard to errors in estimating vortex or source span. It 
is stated in Reference 10 that span-solution within ±0.10B will hold errors 
to an acceptable level. This tolerance is fairly coarse and should not be 
too difficult to attain in practice. An exception, found recently, occurs 
when wall signatures measured in the tunnel corners are employed. These 
locations are significantly more span-sensitive than the central ones. 


Wall Pressure Orifice Location - Peripheral Direction 

Both theoretical and practical considerations arise in selecting the 
peripheral locations for wall pressure orifice rows. Figure 2.^ shows 
theoretical wall pressures, as a function of peripheral location in the 
bound vortex plane and far downstream of a horseshoe vortex in a wind 
tunnel. As expected, center-roof and center-floor locations give the largest 
pressure signals due to lift and so are good candidates, from a theoretical 
standpoint, for upwash interference predictions. The tunnel corner locations, 
5 and 13, are much less sensitive*. While roof locations are usually very 
practical, there may be difficulties with floor orifices. In large tunnels 
there is the obvious problem of foot traffic but in all tunnels powered 
models may involve jets or jet sheets which impinge on the tunnel floor. 

Even if tunnel floor separation is controlled (see Section J*) , jet-impinge- 
ment may compromise the floor pressure signature. 




[l Pressure Orifice Location - Axial Direction 


As indicated^ in Reference 9, Section *4, a test section length of about 
1.5 times tunnel width is desirable to obtain adequate asymptotes to the 
pressure signatures. Orifice spacing should be smallest opposite to the 
model and it's immediate wake and may increase towards the test section 
ends where pressure gradients are less. 


A generous number of orifices should be provided on the floor at and 
ahead of likely jet impingement locations, for monitoring ground boundary 
layer control. in jet- imp i ngement situations, only the forward part of a 
floor-orifice row may be usable for tunnel interference estimation. In 
Other situations a less dense selection from the whole row will be useful. 


Vcrtex-Dnduoed Upwash Effects at the 7unne~ S i dc:ca .. - 

In broad terms, floor and roof orifice rows may be thought of as 
responsible for sensing vortex- i nduced flows and thereby providing data for 
upwash interference corrections. The sidewall orifices are used for 
estimation of blockage corrections. 


Far downstream of the bound vortex, Figure 2.U shows that upwash induced 
by the trailing vortex systems can have a significant effect upon sidewall 
pressures. This could affect the downstream asymptote of the sidewall 
pressure signature and cause an apparent increase in wake blockage. However 
the implicit assumption of Figure 2.^, that the trailing vortex remains 
horizontal, must be reviewed before any legitimate comment can be made 
regarding corrections for such cross flow effects. 


*Mote , 


however, that 13 becomes the proper location for the "sidewall* 1 row 
in ground-effect tests. It is also needed for semi-span model tests. 
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In Reference 13, vortex roll-up calculations are described for wings 
situated in tunnels of various relative sizes and shapes. Though the 
rentral vortex sheet deflected significantly m some cases, the e , 

centers drifted downwards very little. A more extreme exper, mental result 
is presented in Figure 5-2 of Reference 8 concerning flow measurements 
behind a partiai-span jet-flapped model. In this case the 1 ' p V ° r J 6 * . _ 
was horizontal and the flap vortex moved do»n s . go . f leant Jv only . ^ 

C„. On this basis, it appears reasonable to assume that the tra .ng syste 
remains essentially horizontal and to consider correcting sidewall pressure 
signatures for trailing vor tex- i nduced upwash. Since there ,s no cor re - 
ponding source effect on floor/roof increments due to 1'Jt, for ce 1 > Y 
mounted models it is possible to analyze these first and then correct t e 
wall s i qnature * for vortex- i nduced crossflow, prior to setting up the bloc^ag 
^ilySir iio iteration is repaired and the lift and blockage problems rema.n 
essentially uncoupled for unswept configurations. 

2.4 Effects, on Measured Signatures, of Sweep, 

Angle-of-Attack and Model Offset 

At zero angle of attack the addition of sweep to the source and 
vortex lines only affects the shapes of the velocity distributions at t . 
lunnel surface and there are no "cross" effects such as vor tex- , nduced 
apparent blockage or source- i nduced apparent lift. ’ °" * i nter - 

the swept system, these effects appear and must be considered^ To inter 
pret them a relationship must be established between Q/BH and io BH , 
the* respective normalization velocities for source and vortex- . nduced 

effects . 

We may find the ratio of total drag to total lift for a line-source, 
line-vortex system as follows: 


Lift = pU»rb where b is vortex span 
Induced Drag = g- P 


Prof i le Orag = pU.„ Q 


Thus 


’tot ff 


pF^ + p Uoo Q 


p u«rb 


+ A. 

S Ua.b Fb 


2.1 


At (I/O) induced and profile contributions are equal so 


<p 


Thus 


MIN 

20 . . 

fb 


v r 

k Cb 


20 

rb 


<£> 


2.2 


MIN 
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This permits us to interpret the source-s trength , vortex strength relation- 
ship in terms of (L/D)hAX* After some algebra, we obtain 


Q _ , <f> 

r/ BH 2 " (l/d) max 2 - 

For the basic, knee-blown flap model, tested in the 30" x 43" tunnel, 

b/B * 0.465, B/H = 1.433 


so that 


r /BH 


. i*651 x 1 .1972 


TlTdT 


MAX 


.2784 

TlTdT 


MAX 


= .0928 for (L/D) max = 3 
= .0398 for (L/D) max = 7 


2 . 4 


In a typical test case, 
ang 1 e-of-at tack , it was 

tota 1 
and total 
so that 


for the knee blown flap 
found that 

Q/U^BH « 0.0338 
F/UcjBH = 0.5527 
Q/f/BH = 0.0612 


model at Cu = 2 and low 


which is within the range in Equation 2.4. 

Effects of Sweep and Angle-of -Attack on Wall Sig y iatures 

To demonstrate these effects, an example has been selected which is 
based upon the geometry of the swept, knee-blown flap model in the no-tips 
configuration. Effects at the tunnel wall are shown in Figure 2.5* Sweep 
and ang 1 e-of -a t tack effects will be discussed first. 

Figure 2.5(a) shows that adding sweep to the line-source system shifts 
the axial velocity signature downstream. This is expected, since the same 
value of root (X/B) is used. The shift is insensitive to ang I e-of-at tack , 
which is a welcome feature. 

Vortex- i nduced axial velocities at the sidewall (Figure 2.5(b) are 
entirely dependent upon ang 1 e-of-a t tack . For typical relative strength 
values (Equation 2.4) it is apparent that peak 1 ‘c ross M - i nduced velocities 
at high ang 1 e-of-at tack may be comparable with the direct, source- i nduced 
velocities. This probably explains the over-correct ions for blockage 
noted in Reference 9 for swept-wings. 
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Vortex-induced upwash at the sidewall (Figure 2.5(c)) is comparable, 
in norma ized units, with the source induced horizontal velocities 
(Figure 2.5(a)) - as might be expected. It is appropriate to relate the 
upwash to the mainstream velocity: this may be accomplished via the lift 
parameter c l. hl? (i.e. lift coefficient normalized on model span times tunnel 

8 ^ ei9 . t> . 1 ? the P resent case). As a C t_ h b Vd,ue of 2.0, which cor res- 

ponds to incipient tunnel flow breakdown (see Reference 8) the maximum 
value in Figure 2.5(c) of 0.60 represents an upwash equal to about 25% of 
mainstream. When added vectorially to a unit mainstream, an increase of 
only about 5% occurs in the total vector. This would increase somewhat at 
t e higher CL hb values peimissible when ground-blowing is used; correction 
or the effect on blockage is probably desirable at this point. 

Figure 2.5(d) shows that sou roe- i nduced upwash at the sidewall location 
is an order-of-magnitude smaller than the source- i nduced axial velocity 
(Figure 2 5(a). When combined vectorially with the toial axial velocity the 
6 sets o - source- i nduced u pw6s h will be neg 1 i g i b 1 e 

Effects, at the Sideua.ll, of Change to Model Pivot Location 

Curves are included in Figure 2.5 which show the effect of changing 
rom the standard, m i d-sem i -span a-center to one at the wing root. The 
latter was used for swept K8F model tests. At 25-degrees ang 1 e-of-at tack 
this places the entire model approx ima te 1 y 7iS nearer to the tunnel floor’ 

In most of the cases in Figure 2.5, the effects of this change are small.' 

Foi Figure 2.5(d) this is also true because the overall effects are small 
(see above). However the effect on vor tex- i nduced horizontal velocity is 
noticeable and it is apparent that offset effects must be included when 
calculating this correction to the blockage signature. This feature could 
be trouo ! esome because it is ang 1 e-of-a t tack dependent. 


£ fleets on t\Ooj -Mv>iuo-F Loci' Signature 

Figure 2.6(a) shows that the sum of the roof supervelocity and the 
floor counterveiocity, induced by the vortex system, is substantial. Sweep 
reduces the peak velocity differences (u R - up) . It is found that the 
swept vortex curve, at zero ang 1 e-of-a t tack , is essentially unchanged by 
adding 25-degrees of incidence. The pivot location is conseguently 
i mm a t e r i a 1 . 


Source "cross" effects, on the "lifting" ( roof -m i nus- f 1 oor ) signature 
(Figure 2.6(b)) are small when relative vortex/source strength is considered. 
The fact that the forward pivot case produces less "cross" effect is, at 
first sight, surprising. This arises because the tunnel roof and floor 
centerlines are most affected by the central region of the source system 
which remains On the tunnel axis for the forward pivot, but which moves 
towards the roof, with increase in u , for the mid semi-span pivot. 


Ground Effects Testing 

Tor i n-g round -effect testing, either the tunnel floor ("ground") or 
the first ground image may be regarded as part of the model under test, 
he true "center-sidewal I" orifice row is now situated at the foot of the 
t “ nn ®' *' d ® wal1 and > strictly speaking, the blockage sensing orifice row 
should be located here. The roof orifice row remains correctly located but 
in impingement- free cases, the tunnel floor row senses pressures which 
correspon to the wi th-blockage double-tunnel centerline velocity distribution. 
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Though it would be possible set up the necessary computation schemes 
on the above, somewhat idealistic basis, it is more practical to consider 
the in-ground configuration as a beiow-center test when recovering source/ 
sink and vortex strengths from the measured pressure signatures. in the 
second-stage analysis, interference velocities are then calculated at the 
tunnel floor location, rather than at the true tunnel centerline. Both the 
in-tunnel vortex/source arrays and their first ground images are omitted 
when calculating blockage and upwash interference. 

Offset Models 

Some times, the need arises to conduct a "centt^ tunne 1 " test with the 
model displaced vertically from the tunnel centerline. One reason for doing 
this would be to increase ‘ground* clearance so as to reduce the severity 
of impingement problems for powered models. Ground-effects testing would, 
of course, involve below-center models. An orifice row situated at mid- 
sidewall "sees" not only the desired blockage effect associated with (for 
example) an above-center model but also a bound-vcr tex- i nduced counterflow 
which, wrongly interpreted, would appear as a negative solid blockage com- 
ponent. Distortion of the tunnel roof an 1 floor signatures would also occur 
because of offset effects for both vortex and source systems. 

A swept-wing model at ang 1 e-o f-at tack has several similarities to the 
off-center model. The front of the model, situated above-cen te r , has some 
of the properties just described while the tips, below center, yield incre- 
ments of opposite sign and shifted aft. The net effects are illustrated in 
Figures 2.5 and 2.6. 


2.5 Use of Least-Squares Smoothing 


Though the results of the pilot study (Figure 2.3) were encouraging* 
doubts remained about the response of the alternating inverse elements 
(Figure 2.2) to severe data scatter. Figure 2.7 explores this problem. ^ A 
single point on a smooth, 'standard 1 wall pressure s i gna ture > des i gnated ‘A 
in Figure 2.7(a) ? was perturbed upward and downward as indicated at 'B' and 
1 C * „ Though the interference results from case A agreed quite well with 
those derived via the older, iterative method (Figure 2.7(c)) , the conse 
quence of perturbations 1 B ' and ' C 1 were serious (see Figure 2.7(b)). Thi 
provided strong motivation towards a 1 eas t-squa res approach. 

Derivation of Least-Squares Equations 

We define vj* as the velocity induced at the j th observation point 
by the i th singularity and it's image system in the tunnel walls. Due to 
the complete set of N singularities the total velocity induced at the j th 
point is given by 

N 


where oj are the required individual singularity strengths. If the 
corresponding measured velocity Vj differs from the calculated value Tj 
by a residual amount 6j we may write 


or 


6 j = 




H 

■ l 


' J 


o . 

i 


The objective of the least squares approach is to minimize the net 
area between the Vj and the Vj curves as determined at the .V observation 
points. To do this we minimize 
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\ = 1 

0. 
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J 
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, V ' j V kj 

= . 1 v • 
j - 1 J 

V kj 


( 2 . 1 ) 


( 2 . 2 ) 


for 


1 « k < N 


The previous N x N equation set used to obtain an exact match at 
every observation point j is now replaced by an N x N set. N, the number 
of singularities, may be greater than, equal to or (more usually) less 
than N 9 the number of observation points. The case N = N i s not equivalent 
to the "HATCH" procedure described previously because the t heo ret i ca i curve 
is fitted to the experimental data in a 1 eas t -squares sense. On writing 
equation (2.2) in the form 


'V 


to.) = [B. ] 

i k 


we notice that the elements of A., no longer can be identified simply as 
influence coefficients. The elements are no longer simply observed 
velocity increments at k but are now weighted sums of all .V increments. 

Figure 2.8 is the least-squares equivalent cf Figure 2.2, which 
generates an exact match. It is evident that the 1 eas t-squa res process 
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has caused the upper source-sink matrix to become symmetrical about the 
leading diagonal and the largest element is now only 16-times the snal es , 
rather than almost 300-times. However the chequerboard P^s and mmus 
Dattern in the inverse matrix (center table) still remains. The o • 

u!.d to obtain cantor, in. interference direct,, roc «l v.loc.y 
increments, do have a changed structure, however. It may be seen that, 
rather than the previous 'chequerboard' plus-minus structure (r.gure 2.2 
signs now alternate by column. However, the significance of this must be 
appraised via studies of some typical cases. 


Examples of the Least-Squares Approach 

Figure 2.9 repeats the example shown in Figure 2.7, which demonstrated 
the sensitivity of the 'MATCH" approach to data scatter, but applies the 
above, least-squares solution to it. It is evident that the ptevious 
sensitivity to "noise" in the data has been largely eliminated. 

Figure 2.10 is an example of a complex, double-peaked wall pressure 
signature, measured under tunnel flow breakdown conditions with no tunnel 
boundary layer control applied. Though the example is somewhat artificial 
for this reason, it shows that the restriction of the previous, iterative 
method to simple, single-peaked pressure signatures have been removed. 

This flexibility, the data smoothing capability and the reduction of^ 
maxtrix size afforded by the least-squares approach represent a significant 
advance over the previous approaches. 

2.6 Mathematical Summary 


Having reviewed the physics 
previous sections, we are in a 
source and vortex strengths may 
the^e will be set up as direct 


of vortex and source variables in the 
position to set up the equations from which 
be obtained. In the interests of clarity, 
influence rather than least-squares equation 


s . 


Notation 

Subscr i pt 


Summat i ons 


Super scr i pts 


U,V,W 


u,v,w 


C 


P 


i is ilis index for the source or vortex, 
j is the index for the sensing point. 

\ q i 5 the number of source variables, equal to 
the number of wall X- loca t i ons . 

I r is the number of vortex variables, equal to 
the number of roof/floor X-lccations. 

R Roof 
F Floor 

RF Roof value - Floor value 
W S 1 dewa 1 1 

"Direct" influence coefficients, due to unit 
singularity i.e, due to V for rod/ floor 
sensing points and due to Q for wall points. 

"cross'* influence coefficients including both 
axial and norma 1 -to-ma i nst ream effects. 

Measured static pressure coefficient 
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Influence Equattons _ . 

Equating measured and theoretical roof ~m i nus~ f 1 ooi axial velocity 

components, we obtain 

•r 'Q rf 

u. (APT p " - A^T,^ ) - I Ujj r, * l U|J «, 

' t ' \ 

Fig 2.6(a) Fig 2.6(. ^ 


Equating measured and theoretical sidewall pressure coefficients taken 
as the mean of the two sides, we obtain: 


2 'r w 2 Q w 


C pj = u£ - ([ U.. Q.+J f. r.) +( ^ w.. r.) + ( ^ w.j Qj) ( 

t " t t " 1 t 


¥< fi <) 




<4, 

Z SCO 


We note that equations (l) and (2), which will be needed to find 
and Q. are coupled and, because of terms four and five of Equation \2 ) , 
nonlinear. However, we saw previously (Figure 2.5(d) and related dis- 
cussions) that the fifth term is very small. Dropping this term makes 
(2) linear in Q; and permits us to write (l) and (2) as: 


r. = u. . 

' ' j 


Uop (/ 1-C * - / 1 -cJi) - I U. . Q. 



T 1 +T 2 


Q. = U W . 
i 'J 


;u„ (i -Cp ) - ( l 

J i=i 


'v T ,y t 5 




i i 


V u w . r. 

.r, 'J 1 


If the vrl . term is also negligible, as for small span or low lift 
cases. (3) ani J ( U ) may be combined as 



u.O'~ p j - /rn-) 


v'1 - c r 


Solution of Influence Equations . RF 

For the general, large T , case the four sub influence matrices Ujj, 
u?^. U*. and u W . and the upwash matrix W. . are required. The form ot 

equation (5) is J less useful than it appears, not only because it lacks the 
w w. correction but also because data is taken from two distinct populations 
(tire roof/floor and sidwall signatures) which violates an underlying assump- 
tion of least-squares theory. For these reasons an iterative scheme has been 
adopted. This is illustrated in Figure 2.11. c or convenience of layout the 
pressure terms, which are dominant, are shown last in the equations given in 
the figure. 


3.0 TEST MODELS * RIGS AND PROCEDURES 


3.1 General Comments 


Many of the tests which will be described are essentially repeats of 
earlier tests (Reference 8) with augmented wall pressure i nst rumenta t ion 
and, in appropriate cases, floor blowing to suppress tunnel flow breakdown. 
Since detailed descriptions of the models concerned have been given prev- 
iously, particularly in Reference 8, only the ma i n dimensions and details of 
any relevant changes will be presented here* 


The models and rigs 
(subsection 3-2), the un 
3.3) and wind tunnel ins 
30" x 43“ low speed wind 
on the simple wing were 
upwash interference dete 
Selected results are inc 


to be described comprise a simple, semi-span wing 
swept and swept knee~b lown-f 1 ap models (subsection 
trumentat ion. All tests were conducted in the 
tunnel (the “MTF“) at Lockheed-Georg i a . The tests 
conducted as part of an in-house, pilot program on 
rmination by the wall pressure signature method, 
luded in the present report for illustrative purposes. 


3.2 The Simple Wing 

Figure 3*1 shows a f loor-moun ted semi -span wing having a whole-wing 
aspect ratio of 3*0. It has an NACA0012 section and body -of-revo 1 u t i on tips. 
At the quarter-chord location, a 1-inch diameter bar extends downward 
through a clearance bole in the floor and attaches to a 3-component platform 
balance via a turntable which is used tc set ang le-o f -a t tack . The bar may 
be replaced by a cylindrical balance which adds wing normal force, normal 
bending and end load to the lift drag and pitching moment measured by the 
platform balance. There is a clearance of app rox i ma te l y 0.10" between the 
wing root and the tunnel floor. The wing root is immersed in the tunnel 
floor boundary layer which is uncontrolled. Nonetheless, checks between 
data from the present wing and established finite wing theory show minimal 
performance degradation due to wing root effect. 

The photograph of Figure 3*1 was taken through a new, laser velocimeter 
window which now comprises the back wall of the 30“ x 43" test section. Part 
of the laser velocimeter may be seen at the right. 

3*3 The Unswept and Swept Knee -B 1 own- F l ap Models 


Figure 3*2 was also taken through the new back window/wall of the test 
section. Though the swept knee-blown- f lap model is the object of the photo- 
graph, a good view of the sting, model air supply, tunnel wall pressure orifice 
strips and the floor blowing slot are also obtained. Though the sting appears 
quite massive in this view, it should be noted that It is only about 2-inches 
wide. Most of it disappears into the floor at high ang le-of-a t tack , as shown 
in Figure 3*3* 


Figures 3*^ and 3*5 show the principal dimensions of the unswept and 
25“degrees -swept knee blown flap models. For both models the tips and the 
slats are removable. The flaps are integral with the model, however and 
have upper surface angles of 76- and 60-degrees to the wing reference line 
respectively for the unswept * and swept -w i ng models. Further dimensional 
and sectional details are given in TABLE I. 
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3.4 Floor Boundary Layer Control 


The boundary layer control rig used for ground blow-ng in previous tests 
(References 3* 6 and 8) was modified for the present test series by providing 
the capability to control three spanwise slot segments independently. Sep- 
arate controls were provided for a central 8-inch span slot and two 6-inch 
segments to each side of this, for a 20-inch total, equal to the powered span 
as recommended in Reference 8. (Previous tests employed a 30-inch span slot). 
The change in supply arrangements made it necessary to revise the blowing slot 
detail to the form indicated in Figure 3-6. The slots had been situated above 
the middle of each plenum in earlier tests. Spacers were used at regular span- 
wise intervals to maintain the 0.067~inch slot height. More were required than 
previously because of a change from stainless steel to aluminum plenum covers. 

Slot calibration procedures were as documented in Reference 3- As before, 
blowing rate was monitored using plenum static pressure taps. 


3-5 Wall Pressure I nst rumen tat ion 

Figure 3-7 shows details of wall pressure orifice locations used for 
tests on the knee-blown flap models. it should be noted that, for these tests, 
rows 3 and 5 were located on the upper and lower side walls and not on the 
roof and floor as shown in Figure 3-2. The orifice strips were moved after 
completion of the main tests to accommoda te the laser velocimeter window. 

Previous i ns t rumen ta t ion comprised the sidewall orifice rows, 2 and 4 and 
the floor rows, 7 and 8. The latter rows were augmented for the present tests 
to give better resolution for identifying the ground vortex and hence flow 
breakdown. Rows 1, 3 and 5, in the tunnel corners, are new. Rows 1, 3, 5, 6 
and the aft parts of 7 and 8 were made from aluminum strips, as may be seen in 
Figure 3-2. This, newer arrangement is preferable to orifices installed 
directly in the tunnel walls. General comments about pressure orifices, their 
location and their use may be found in Section 4 of Reference 9- 


3.6 Tunnel Speed Control 

The desirability of running at "cor rec ted-q H during powered model tests 
is well known. In previous tests in the present series (References 6 and 8) 
this was achieved by sensing wall pressures upstream and downstream of the 
model at suitable locations and using a voltage divider network (Figure 3.5 
of Reference 6) to interpolate for an effective pressure at the model location. 
Though this approach was quite successful, the fact that it relies upon only 
two pressures, rather than a whole pressure signature* is an obvious weakness. 

A specific shortcoming is that solid or separation bubb 1 e - i nduced blockage is 
likely to be underestimated. 

For the present tests, the matrix method for blockage was available in 
time to permit on-line, whole-signature analysis to be used for speed control. 

A combined inverse and centerline interference matrix (similar to Figure 2.7, 
lower part) was applied to supervelocity data derived from the sidewall orifice 
rows. Tunnel 'q* and thereby C u was determined at the model using the on-line 
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data reduction program and the tunnel speed control was adjusted until the 
desired C u was obtained. At the time of testing, no swept-bound vortex cap- 
ability had been developed, so a straight wing matrix was used for the swept 
wing tests. 
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4.0 U^L OF TUNNEL-FLOOR BLC TO SUPPRESS FLOW BREAKDOWN 


4.1 Effects of Tunnel Blockage and Flow Breakdown. 

The major problems confronting the test engineer in a powered model test 
have been, in order of decreasing importance: the difficulty in running “whole" 

C u ‘s, the related difficulty in correcting forces for blockage effects on 1 q 1 , 
the difficulty in recognizing when flow breakdown effects have become excessive, 
the impossibility of correcting for them and, finally the problem of angle-of- 
attack correction with curved, powered wakes present. 

It is believed that the studies described below represent the first 
successful attempt to solve the overall problem and identify the specific con- 
tributions, to model forces, attributable to the various effects mentioned 
above. The general approach will be to start with uncorrected lift data for 
the unswept knee-blown flap model, at high C u , and illustrate the ef f ects of 
first correcting for blockage and then applying floor blc to suppress tunnel 
flow breakdown. 

Blockage and Angle -of- At tack Corrections 

Figure 4.1 shows *■ a curves measured using on-line blockage corrections, 
as described in Section 3, at 'whole' C u values of 4.0 and 10.0. For comparison, 
"free air" curves are included (broken lines) which represent data measured in 
the 7 ' x I0‘ tunnel at NASA-Ames. The crosses in Figure 4.1 show C|_-values 
which employ nominal tunnel -q and uncorrected v values. Since corrected - C ;> 
is held constant, uncorrected values vary with n and are greater than the 
set values. 

The circles represent data corrected for blockage, by the matrix method 
and for ang le-of-at tack , by the Williams and Butler method (see References 11 
and 6 - section 5)* We shall see in Section 5 that use of pressure signatures 

to determine ang 1 e-of -a t tack correction procedures almost identical results 
in many cases. Chained lines in Figure 4.1 connect cor respond i ng uncorrected 
(crosses) and corrected data points. It is evident that although corrections' 
are reasonably successful at lower ang 1 es-of-at tack and C. p -values, signif- 
icant errors remain at high a's. 

Use of Tunnel -Floor Blowing 

The first tests on the unswept knee-blown flap model in the present 
series were Jsed to develop ground-blowing strategy, recognizing that, in 
distinction to previous tests, the objective is to remove the ground vortex 
entirely, if possible. For the previous, g round -e f fee t tests the objective 
was to establish a zero-skin friction condition at the ground. 

Several candidate criteria were considered for determining the tunnel- 
floor BLC setting. However, it rapidly became apparent that the best procedure, 
was to eliminate entirely the negative pressures upstream of the jet impinge- 
ment, as illustrated in Figure 4.2. Line printer symbol plots were made 
routinely of the center-floor static pressure signature and blowing was in- 
creased until the suction peak disappeared. No attempt was made to prevent 
jet impingement. At this point, no force correlations had been made and no 
flow field measurements had been attempted. 
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Figure 4.3 shows that the use of f 1 oor -bl owi ng to suppress flow break- 
down was remarkably successful in removing the residual errors in the Figure 
4.1 lift curves. The errors in the previous blockage and incidence corrected 
data (circles), were virtually eliminated when floor-blowing was applied 
(triangles). Only for the last two points at C u * 10 was f 1 oor-b 1 owi ng not 
fully effective: for these the limit of blowing capability evidently had been 
reached . 

Figure 4.3 also demonstrates the significance of the distinction between 
ground boundary conditions appropriate to ground-effects as opposed to center- 
tunnel testing. The moving-ground points (pluses in Figure 4.3) give the 
correct result for a ground-ef fee ts case. It is evident that, for this case, 
a floor vortex should be oresent , rolling just above the moving ground. Because 
of this, some lift degradation would occur, relative to the cor respond i ng 
free air case* in flight near the ground. 

Magnitudes of correction and ground blotting quantities 

Figure 4.4 shows typical blockage corrections, ang 1 e-of -a t tack corrections 
and ground-blowing C ;j values as a function of ang I e-of -a t tack at typical model 
C u values. 

In the most extreme case, the tunnel-q setting was only 65% of the q 
experienced by the model. Ang I e-of-at tack corrections appear to be less 
sensitive to C u and peak at about 4-degrees. Some scatter is evident in the 
ground-blowing C u settings but the general trends are clear. Though values 
of the order of 0.6, for the C ^ - 10 case, seem high they are a small fraction 
of the corresponding model blowing momentum coefficients. The blowing 
pressure ratio scale, to the right of the ground C M plot in Figure 4.4 does 
net apply to the C u = 10 case because this was obtained via a reduction in 
tunnel-q at constant model mass flow. 


4.2 Flow Measurements Using the Laser Velocimeter 

At the end of the planned test series on the knee-blown flap models, 
tunnel modifications were made to install the .arge window in the back wall 
of the test section* The laser velocimeter was then installed, in preparation 

for another test. The opportunity was taken to investigate the flow breakdown 

phenomena just described by making LV flow traverses near to the model center 
plane. To reduce "shadowing’', the s t ra i ght -wl nged model was reinstalled. 

Figures 4.5 and 4.6 provide vivid evidence of flow breakdown, ahead of 
the model, at extreme model C u ‘s and angles of attack. These are fixed- 
floor cases with no blc applied. It will be noted that the incident flow 
angles, just ahead of the model are quite low in Figure 4.5 and 4.6. Figures 

4.7 and 4.8 show the same model conditions with blowing applied at the tunnel 

floor; The floor vortex has been pushed back almost to the impingement point 
in both cases and it's size has been reduced markedly. Just ahead of the 
model, the incident flow angles are much greater and the flow vectors are 
longer. These changes are consistent with the lift increases observed when 
floor blowing was applied. 

The data of Figures 4.5 through 4.8 confirm the choice of the criterion, 
discussed previously, of increasing floor blowing until suctions below the 
floor vortex vanish. 
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4.3 Interpretation of Wall Pressure Signature^ 


Wall supervelocity detc., derived from pressure signatures for 
previous = 6.0, a = 28° case, are shown in Figures 4 .9 and 4.10. 
floor blowing (solid points) the floor vortex peak is readily ident 
rows 7 and 5 and may contribute to the row 2 and 4 (i.e. sidewall) 
However, it appears that row 3 is not affected: its peak is too far 

to be vortex-related. 


the 

Wi th no 
i f i ab I e in 
peaks . 
forward 


slot (+ symbols), the suction 
row 7 (by de f inition), but is 


On applying blowing through the 20-inch 
peak disappears entirely at the center floor 
not entirely removed at row 5, the lower sidewall, where there is no blowing 
to suppress it. It is also very likely that the second peak in rows 2 and 4 
also marks the path of the ’’floor" vortex in the 20" blown-floor case. 


Data with the slot-width reduced to 8 inches (triangles) shows that this 
is less effective than the standard, 20-inch slot. This is confirmed by force 
measurement s . 


E^'eets on Tunnc" Corrections 

It is disturbing that, under high-C p , high-a conditions, the main 
suction oeak measured at the sidewails (rows 2 and b) may include a signifi- 
cant component caused by the passage of the floor vortex across the tunnel 
sidewall orifice row. However, the application of floor blowing shifts the 
vortex aft leaving what is probably the correct, solid/bubble-blockage induced 
first peak. Nonetheless, the wall pressure signature input to the tunnel 
blockage correction program still includes a second peak which is directly 
induced by a vortex, rather than being a true reflection of tunnel blockage. 


Figure li.11, taken from Reference 10, explores the effect of a dominant 
second peak upon tunnel centerline blockage interference. The experimental 
case (circles) is compared with an idealized case (triangles) with the 
second peak removed. As might be expected, there are significant changes in 
interference opposite to the second peak itself. However, the effect of the 
peak on the interference at the model location is surprisingly small. 


Another effect of floor blowing, indicated to some degree in the row 3 
data of Figure 4.10, is a general reduction in blockage interference. This 
is especially noticeable on applying BLC to the C u = 10 case illustrated in 
Figure 4.12. In this case, it is speculated that, with no BLC applied, the 
separation streamline, from ahead of the floor vortex, rises to perhaps halt 
the model altitude at its crest. Though the pressure signature blockage 
prediction method probably responds to this with appropriately located cor- 
rections of the proper sign, the total flow is far too distorted for any 
such corrections to be taken seriously. It is obviously better to get r.d . 
of the floor vortex, by applying blowing blc, than to try to correct for it s 
effects. 


Addition: Floor-Vortex DatJ 

Figures 4.13 and 4.14 show the development of floor centerline 
pressure distributions, with angl e-of-at tack , at “ 2.0 and C„ » 4.0 
respectively, with no floor blc applied. Figure 4.15 summarizes the data 
for Cy * 4.0 in terms of vortex and impingement location. Correspond i no 
pressure data, at both C w 's are presented in Figure 4.l6. 
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The impingement point moves forward, as expected, with angle-of-attack 
(Figure 4.15). The maximum suction point remains an almost-constant distance 
ahead of impingement which suggests that vortex size is not very dependent 
upon model angle-of-attack. The first positive pressure peak gives a general 
indication of the location of the ground separation point, however the peaks 
are not well defined at the high angles of attack. 

The development of peak pressures is shown in Figure 4 . 1 6 . At Cu = 2 
this plot gives a good definition of the angle-of-attack for the onset of floor 
separation, i. e. where the ''ortex and impingement curves diverge from the 
single, first positive peak, line. 

Application of floor blowing eliminates the vortex suction lines in 
Figure 4.16. However there is very little change in the impingement pressure 
curves . 
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5.0 


FORCE AND MOMENT CORRELATIONS 


5-1 Checkout for a Simple Wing 

Before embarking upon an investigation of powered '-mode ! corrections, it 
appeared desirable to test the new, wall pressure s i gnature-based angle-of- 
attack correction procedure on a simple model. Appropriate tunnel pressure 
and model force data were therefore obtained in tests on the wing shown in 
Figure 3.1. 

Figure 5.1 compares Cj^ - a curves corrected by the classical, Glauert 
method (-^-symbols) and by the new wall pressure signature method (circles). 
Wall pressure s i gna ture - de r i ved blockage corrections, which were small, 
we re applied in bo th cases. It is evident that the new me t hod p rov ides angle- 
of-attack correction estimates slightly smaller than those determined via the 
‘Glauert 1 approach. However the generally good agreement gave confidence 
that the new method works properly. 

5.1 Selection of Singularity Geometry and Iteration Procedure 

Geometry of Vortex and Source Elements 

The effects of sweep and ang 1 e-of -a t tack on tunnel influence coefficients 
were discussed in some detail in Section 2 from a theoretical standpoint. 
Figure 5.2 shows results from a practical application to the swept-wing, 
knee-blown, jet-flap model in a test at C u = 2. Influence matrices corres- 
ponding to two different element geometries were used in Figure 5-2(a) to 
correct measured data (chain lines) for comparison with "free air" data 
(circles) measured in the NASA-Ames 7 1 x 1 0 ' wind tunnel. 

The broken lines show corrections based upon the "correct" swept element 
geometry set at 15“degrees ang 1 e-of -a t tack for both bound vortices and 
sources. For the full lines, simple, unswept elements were used. the 
influence of geometry is clearly very minor for this model and tunnel com- 
b i na t i on . 

Effect of "cross’* terms 

Figure 5.2(a) displays over-correction of the lift curve. It has al- 
ready been mentioned that previous blockage over-cor rec 1 1 on may have been 
a consequence of neglecting the effects of vor tex- i nduced upwash "cross" 
effects on the sidewall signature. The broken lines in Figure 5.2(b) show 
the results of a full iteration, as outlined in Figure 2.11, applied to the 
previous example. The differences between the broken lines and the full 
lines are the effects of "cross" terms. As anticipated, the over-cor rect i on 
of the lift curve has been almost eliminated. 

Fur the i examination of the "cross" terms revealed that w-squared term - 
i.e. wake upwash at the tunnel wall centerline is by far the largest. It 
is also found that the effects of these terms become excessive beyond C v _ - 2 
(see Figure 5.2(c)). TH s suggests that the horizontal trailing vortex 
model starts to fail because its 1 geometry is fixed. It may be shown that 
incremental C^-cor rect t ons due to "he effect of w-sqi.ared on blockage 
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proportional to C^-cubed for a fixed-geometry wake. The system is very 
sensitive, at high to small changes in wake location, 

5.3 Analysis of Ang 1 e-of -At tack Corrections 

The corrections for the knee-b 1 own- f 1 ap models are dominated by blockage 
effects and the sensitivity to errors in ang 1 e-o f -a t tack corrections is quite 
small when plotted in conventional lift curve and drag polar form. me 
present angl e-of -at tack corrections will therefore be assessed in comparison 
with other predictions. 

The wal 1 pressure signature method provides a continuous distribution 
of A a along the tunnel axis and an effective model position must be selected 
which cha rac te r i zes its aerodynamics. In the present case, a fixed location 
at x = 0 has been selected for both unswept and swept wings, recognizing 
that other locations - such as varying load centers derived from and C|_“ 
could be considered. 

For the swept wing, the choice of the correction location at the root 
quarter-chord could be questioned. However an aft shift in C.P. location 
on adding sweep did not occur because a lower flap angle was also introduced, 
during design, to improve the drag polar. In fact, the measured swept-wing 
C.P. lay slightly forward of that for the straight wing in most cases. Both 
lay between the quarter and three-quarter root chord locations and moved for- 
ward or aft depending upon the balance between wing and flap lift. 

Figure 5-3 shows ang 1 e-of -a t tack corrections, 4 \ for the basic swept 
wing as a function of blockage-corrected C^. The three parts correspond to 
a) unpowered or low-C^ (i.e. BLC) conditions, b) moderate C^'s with no 
floor impingement and c) cases with floor impingement, with floor blowing 
used. In all cases full-length roof and floor pressure signatures were 
employed, recognizing that errors arise from impingement regions. 

Figure 5.3(a) shows that, as for the simple wing, wa 1 l -s i gnat ure 
derived ang 1 e-of-a t tack corrections are slightly lower than the classical 
G I auert method but increase, per C^_, during and after stall. The hook- 
shaped A® curves occur because the wing cen ter-of-pres sure moves back less 
rapidly with C|_ after stall than it had moved forward prior to stall. It 
is known that the flap separates before the leading edge does at zero and 
1 ow Cp va 1 ues . 

At moderate values, with no floor impingement. Figure 5*3(b) shows 
smaller ang 1 e-of -at tack corrections than both the Glauert (straight-wing) 
and the Wi 1 1 i ams-Bu t 1 e r estimates, even though the latter include a C* - 
related attenuation factor. Howt /er the increase with Cl is more rapid for 
the signature-based estimates than for the others- Though the C.P. does 
move forward with ang 1 e-of -a t tack in both the cases shown^the streamwise 
ang 1 e-of-a t tack gradients are insufficient for this to be the full explan- 
ation. Changing flow geometry may also be partly responsible. This is 
almost certainly true in the impingement cases, at C VJ = *4.0 and 6.0, shown 
in Figure 5.3(c). Here, the trends are generally similar to those of the 
previous figure but the levels are in better agreement with the other 
es t i mates . 


Angle-of-attack corrections for the unswept KBF model (not shown) are 
generally greater than for the swept geometry as should be expected, and lie 
above the Glauert values increasingly up to C u ® 2 . Above this, difficulties 
in signature analysis obscured the trends. 

Comments 

In the non- impinging cases described above, the Aa estimates by various 
methods are generally within a spread of about one-degree. Within this 
range, there is no experimental basis for saying which resulc is correct. 

Further refinement would probably require investigations of surface pressures 
- particularly leading edge suction peaks, in large and small tunnels. As 
mentioned previously, the consequences of these differences to the present 
force and moment data are ^ ^ior importance. 

For impinging cases, the signature method indicates quite large Aot 
values at high compared with the simpler theories. However, the theore- 

tical model used n these cases is clearly inadequate because it fails to 
recognize impir ^ Improvement to the correction orocedure is also 
required in these cases with reqard to the roof/floor pa. t of signature analysis. 

5. A Force and Moment Correlations for Unswept- and Swept-winged knee-blown 

jet- flap models. 

Figures 5.A(a) through 5.7(c) show "free air" (broken lines) and 
corrected small tunnel force and moment data (points) for the four model 
con f i gu ra t i ons tested. In analyzing the data, 'whole 1 floor signatures were 
used in all cases except the straight wing at C u 's of A and above, which 
failed to converge using this procedure. In these cases, the roof signature 
only was used for anal e-of -at tack correction, after removing blockage effects 
and doubling the roof per turba t ions . For the same reason, computation was 
stopped after the first pass for all configurations when C u was A or greater. 

All uncoupled solutions (i.e. independent ang 1 e-of -a t tack and blockage 
solutions) are designated by an asterisk in the C u table. 

NASA CR 152, 2A1 (Ref 9) documents the first attempt to appl\ the wall 
pressure signature method to the present configurations. Relative to the 
earlier, 'q-pot 1 corrections of CR 152,032 (Ref 8) the Ref 9 pressure sign- 
ature results were disappointing because the correlation with "free air" 
data was significantly worse for high-C u cases. This occurred largely 
because of flow breakdown itself but also to some degree because the signa- 
ture analysis of the Ref 9 iterative method can respond adequately only to 
classical, single-peaked pressure signatures. 

Straight-Wing y With and Without Slats 

The previous over-correction tendency of the iterative, Reference 9 
method has been largeriy overcome as a result of the better flex.bility 
of the present method. The chance for success is increased further by 
ground blowing, as was seen in Section A. Figures 5-^ and 5.5 show that 
the present method improves upon both the Ref 9 and the Ref 8 approaches. 

The latter had a tendency to over-correct at high C u and under -cor rect at 
low C u . The overall agreement is now within the limits of experimental 
error . 
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Selected unblown ground data have been added to Figure 5.4 (flagged 
points) to supplement Figure 4.3 which is based upon interim blockage 
corrections and includes Wi 1 1 i ams-But ler ang 1 e-of -a t tack corrections. 

The crosses in Figure 5-4, at C p = 1 , correspond to an "overblown" 
g round-b lowing case in which the blowing was set as for C u = 6, a = 30. 

The C u = 1 case does not include impingement, so the results show that 
blowing maybe left operative at a "set and forget", worst-case level without 
signiricant change to other data. Any q-changes, due to excessive floor 
blowing, are accommodated automatically via the wall pressure signature 
blockage correction procedure. J 

Swept Wing With and Without Tips Fitted 

Relative to previous methods for correction the present swept wing 
corrections, Figures 5.6 and 5.7, show definite improvements in the drag 
polar correlations. However, lift curve slope still appears to have been 
over-corrected at or above = 4, particularly for the with-tips case 
(Figure 5.7(a)). Pitching moments are less well corrected at C u = 6 and 10 
with the tips added (Figure 5.7(c)) but continue to agree well for the 
basic swept case (Figure 5.6(c)). 

With the above relatively few qualifications it appears that the 
differences between the corrected and the free air data are not only within 
the experimental error band but have reached the point where possible cor- 
rections to the large-tunnel data should be reviewed. Rough calculations 
indicate that a large-tunnel C[_-value of 10 would be reduced by approxi- 
mately 0.2 on correcting for blockage effects. This is of the same magni- 
tude as the anticipated experimental error. 
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6.0 DISCUSSION 
6 * I Aerociynam \ cs 


Overview 

The methods described above bring to the wall pressure signature method 
new, more powerful and more comprehensive capab i I i t i es . These include an 
oroer-of -magn i tude reduction in run time and ang 1 e-of -a t tack correction cap a~ 
bility. During development, an effort has been made to make an effective 
trade between flexibility and ease of use. Some typical simplifications 
which have been included are the use of ‘’whole 11 tunnel floor pressure signatures 
knowing that they contain impingement spikes, and the use of a constant effec- 
t i ve model location, rather than one which responds to known changes in C. P. 
location. Despite these self-imposed restrictions the present methods have 
achieved good successes. 

There were major questions, at the start of the present work, concerning 
the sensitivity of the wall pressure signature method to model geometry, 
particularly to sweep and ang 1 e-of-a t tack effects. Strong sensitivity would 
have made the method much less useful. In most respects i nsens i t i v i ty has 
been found not only to sweep and ang 1 e-of -a t tack but a 1 so to sT ng u 1 a r i t y 
spacing. Sensitivity has been found, however, to vortex wake location under 
high lift conditions. This will be discussed below. 

Cases with jet- 1 mp i ngemen t and floor-blowing were the subject of an 
extension to the work planned originally. The use of floor tangential blowing 
and wall pressure signature based blockage corrections is a prerequisite to 
several of the discussions of impingement effects at high-C^ which appear 
elsewhere in this report. However, theoretical modeling for cases with 
impingement is currently much less advanced than for cases without it. It 
is anticipated that j e t- i n-c rossf ! ow experience (see Part II) will help sig- 
nificantly when improving impinged-jet flow models. Further discussion of 
impingement modeling and its problems will follow that for non- imp i ng i ng cases, 
below. 

Wake Modeti>%g for Non-Irnpinging Cases 

AGARD Report 692 is a country-by-country review of 'wind tunnel correction 
methods for high ang 1 e-of -a t tack models. Repeated reference is made to the 
fact that tunne 1 - i nduced wake distortion must be considered, even for unpowered 
models, during the correction process. This appears to contradict the assertion, 
in sub-section 2.3, that wake vortex movement is not significant for the 
present KBF tests. While the present very high~C^ test results support the 
AGARD 692 assertions, an apparent paradox remains in the lower ranges, in- 
cluding most of the region of practical interest. 

A vortex fair shed to the 50-percent semi-width positions in a rectangular 
wind tunnel, at mid-height, possesses the special property of being in equili- 
brium in the cross-plane. This may be confirmed by considering image vortices 
which give cancelling induced velocities at say, the right-hand trailer (see 
Figure 0.1). Members of a vortex pa i r shed near to these special points in tne 
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tunnel orbit them at a rate determined by vortex strength. This appears to 
be the situation for the basic KBF models, which span 46.5% of the width of 
the 'MTF' tunnel. This choice of model size may explain the paradox men- 
tioned above. 

The above special result for 50% semi-span models is not new, but it 
may have new significance for sizing and positioning models in tests at very 
high lift. For example, if a powered model's span must exceed 50% of the 
tunnel width, cons i dera t i on of near-wake distortion could be used to locate 
an optimum, above-center location in the wind tunnel. This would also 
relieve the impingement problem. Further work in this area appears worth- 
wh i 1 e . 

Wake Modeling for Impinging Cases 

The most significant property of an impinging flow of the present type 
is probably that the circulation and span of the downstream vortex wake no 
longer define the total lift on the model. This is because the existence 
of the floor stagnation point permits vortex lines to link to the tunnel 
floor. Mathematically, the connectedness of the region is changed by floor 
stagnation and closed circuits can no longer be drawn which define the model's 
bound vorticity. It is partly for this reason that an obvious step, of 
linking total vortex strength to model lift, was not incorporated as part 
of the ang 1 e-of -a t tack correction procedure in the present work. 

The first problem encountered in setting up an impingement model is 
to determine how much of the bound vorticity trails downstream and how much 
joins to the tunnel floor via the jet. In unb 1 own-f 1 oor cases, the standing 
floor separation vortex (Figure 4.5 etc.) tends to confuse the issue. Some 
tentative trials have been made using tunnel surface pressures, in blown 
floor cases. Tests with several combinations of mode 1 - to-f 1 oor and trailing 
vortex systems gave disappointing results. Tests at C u = 6 and 10, and 
usually 4.0, (Figures 5.4 through 5-7), include impingment, but were cor- 
rected using flow models which ignore it entirely, except via the fact that 
tunnel surface boundary conditions are satisfied. The success of the 
corrections is difficult to explain. While some necessary conditions for 
this success certainly can be identified in the present studies (for 
example, the use of floor blowing and s i gna t u re- t ype methods) additional 
conditions are needed to complete a sufficient set. These are difficult 
to identify: the topic requires further work. 

6.2 Signature Analysis 


Inverse Methods 

Recognizing that the first stage of signature analysis comprises the 
solution of a three-dimensional inverse problem, the methods in Section 2 
were reviewed by a researcher in inverse wing design. It was found that 
conventional inverse techniques could be applied to the present problem. A 
paneled shape cor respond i ng to the model and its wake might be found using 
tunnel wall pressures as the objective function, leading eventually to inter- 
ference velocities. Further review reveals, however, that this approach 
would neither be sufficiently compact, nor sufficiently fast for practical 
use in routine wind tunnel correction work. 
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Flou) Model Gconctry 

Though the insensitivity of the overall method to most model details, 
mentioned above, is interesting ae rodynami ca 11 y , its predominant importance 
lies in the simplifications it affords when the new methods are applied. If 
significant sensitivity to sweep and/or ang 1 e~of -a t tack ^in particular) had 
been found, individual influence matrices might have been needed for (at 
worst) every data point. As it has transpired, relatively few matrices will 
be needed for any particular test. 

While introduction of ,, oross ,, ‘*term capability (Sections 2 and 5) lead 
to the above result, it also revealed sensitivity to wake location as dis- 
cussed above. Within the present framework, the effect has been to limit 
quite severely the use of “cross* 1 - te rms to improve the results. New ex- 
perimental and/or theoretical techniques are needed, to locate the vortex wake, 
before the capabilities of the present methods can be exploited fully. 

Da ta r Condi tio>\ i no ' 

As for the previous, iterative method, the main task of data pre- 
conditioning is to subtract empty tunnel wall super ve 1 oc i t i es from corres- 
ponding mode 1 -present data (see NASA CR152, 24l, Section 4). A subsequent 
conditioning task, in the preent case, will be the removal of jet- in-crossf low 
induced velocity components when appropriate. The last conditioning stage, 
which concerns data smoothing, is embedded in the signature analysis itself 
both for the previous, non-iterative and for the present method. Though 
the latter employs a least-squares procedure for signature fitting, recent 
experience has shown this to be insufficient to prevent a blocked pressure 
orifice at the ‘wrong 1 location from spoiling otherwise good data. A bad- 
point rejection filter, similar to that used in the earlier algorithm, is 
needed . 

The impingement-case f 1 oor- s i gna t ure is a heavy candidate for data con- 
ditioning. Currently, it is either accepted in full* at user option, or it 
is rejected in favor of ’doub!ed-up‘ tunnel roof data. This is not always a 
good alternative. However, there are a number of unanswered questions con- 
cerning how an impingement ‘spike’* for example, should be treated (e.g. : 

Is it theoretically correct to f a i r it out?). The answers to such questions 
should become more apparent as impingement modeling becomes better understood 
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7-0 CONCLUSIONS 


Recent advances in wall pressure signature methods are described and 
used to estimate ang 1 e-of -a t tack and blockage constraint effects for several 
powered models in low-speed tunnel tests. Tunnel floor BLC was employed at 
high C u , to control flow breakdown. The combined techniques permitted 
successful testing well beyond usually accepted limits. 

Use of Tunnel-Floor BLC to Control Flou Breakdown 

Control of tunnel flow breakdown was accomplished using tangential 
blowing, along the tunnel floor, from a point just ahead of the model (Figure 
3.6). Floor pressures were monitored to determine blowing settings. Sub- 
sequent f 1 ow me a su remen ts i with an L V a nemome ter, s howed t ha t the f 1 00 r BLC 
had destroyed the vortex ahead of jet impingement. 


Other observations include: 

(1) Elimination of the floor vortex resulted in a large increase in 
upwash at the model location. 

(2) Lift loss relative to 'free-air 1 conditions at high C u was eliminated. 

(3) There was a significant reduction in tunnel blockage when floor 
BLC was used. 

(k) Floor pressures may be used to monitor vortex destruction: floor 

blowing is increased until the suction peak ahead of impingement 
is eliminated. 

(5) Overblowing is not harmful. An entire test may be performed, 
without detriment, with blowing set for "worst -case' 1 conditions. 

(6) The span of the blowing slot must be no less than the powered span 
of the node I . 

(7) The BLC needed in the present tests significantly exceeded that 
which would be provided by a moving ground matched to tunnel speed. 


Use cf Wall Fressure Bioraiures for Anole-of-At took Correetioi 

Tunnel roof and floor pressures were used to determine the strengths of 
horseshoe vortices* used to represent model lift effects, and thence anqle- 
of-attack corrections. The technique was very successful for a simple wing 
but the corrections for powered models were less easy to interpret because 
in many cases strong blockage effects and floor impingement were also present. 
Some specific observations: 


(8) Aa estimates for unpowered cases and for low-range powered cases 
were generally slightly lower than the classical 'Glauert* 

pred 1 ct i ons . 

(9) in the low- to medium-Cu range, A.t values were comparable with 
Will tams/Buf ler estimates (Ref. 11) at low ang 1 e-of -a t tack but 
increased more rapidly with a . 
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(10) At nigh-C^, Aa values determined from roof/floor pressures were 
generally high. Other tendencies were as just noted. 

It is not possible to judge, from the present experiments, which of the 
Aa estimates was ‘correct'. it is possible that the ‘changing-6' effect is 
related to changes in jet-sheet geometry as ang le~of -a t tack is increased. 

Combined Blockage and A nj I c -of -A ttack Correct ten Program 

New developments in the wall pressure signature method include: 

O Angle of attack correction capability (see above). 

o Use of fixed geometry, mu 1 1 i - s i ngu 1 ar i ty solutions which replace the 
previous, iterative mov i ng- s i ngu 1 ar i ty procedure (Ref. 9). 

o Application of a 1 eas t - squ* r es approach which gives both smoothed fits 
to experimental data and reduced matrix size. 

o Use of a generalized singularity routine tc generate influence 
coefficients for swept geometries at ang 1 e-of-at tack , including 
"cross'* effects for non-planar cases. The latter are also applicable 
to offset mode 1 s . 

o Generation of a combined blockage and ang 1 e-of -a t tack algorithm cap- 
able of handling non-planar "cross" effects (see Figure 2.1l). 

The following facts have emerged: 

(11) The matrix method is almost an orde r -of -magn i tude faster than 
the previous iterative method when applied to a given problem. 

(12) The least- squares approach works well for smoothing ‘local noise' 
but an additional po i n t - re j ec t i on scheme is required for “rogue" 
points (blocked orifices, electrical 'spikes' etc). This has not 
been implemented. 

(13) The least -squares approach cannot and should not be used when 

combining blockage and ang 1 e-of -a t tack solutions: iteration 

between these is effective. 

(1*0 The use of swept singularities, at ang 1 e -of -a t tack produced little 
change relative to corresponding s t ra i ght -w i ng results in the 
present applications (Figure 5.2(a)). 

(15) For the cases investigated, the only significant coupling between 
lift and blockage solutions was via t ra i 1 i ng- vor tex- i nduced upwash 
on sidewall blockage signatures. 

( 16 ) The aoove coupling is significant at high Cj : incremental 
(^ ' s , due to wall upwash, are proportional to -cubed. 

The coupling term is ve» / sensitive to wake geometry. This 
was a limiting factor in the present application. 
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fications should be noted: 

(18) Fully coupled solutions gave improved results, particularly in 
reducing overcorrection for blockage, only for C u less than or 

equal to 2. 

09) Above C.. = 2. the inclusion of the sidewall upwash term reduced 
the blockage correction too much. It is suspecte t at e 
assumption of an undeflected vortex wake may be responsib e fo 
this at mid-range C u 's. Improved imping i ng-jet flow models 
required for high-C^ cases. 

/ ^ \ tendency to overcorrect )ift curves at high C ; , noted 

(20> p«v ousw for toe coof ipurat loo with tips fitted (F,,urP 5 .7(,) . 
has been reduced but not eliminated by the present methods. Drag 
polar and pitching moment correlations are quite good. 
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TABLE I 


MODEL DIMENSIONS 


Fuse 1 age : 

I eng th 

maximum width 
maximum height 
maximum cross- sec t i on 
equivalent diameter 
nose location 
fineness ratio 

balance centerline location: 
water 1 i ne 
butt line 
reference point: 

fuselage station 
water li ne- 
butt line 

Straight Wing : 
sweep 

quarter chord MAC location: 
fuselage station 
water line 
butt 1 i ne 


31. 

55 

cm 

(12.42 

i n) 

4. 

46 

cm 

(1.76 

in) 

7. 

76 

cm 

(3-06 

in) 

30. 

30 

cm- 

(4.70 

in 2 ) 

6. 

21 

cm 

(2.44 

in) 

0. 

00 

cm 

(FS 0.00 

in) 

5 

08 


5-08 


40. 

,6b 

cm 

(16.00 

in) 

0. 

.00 

cm 

0.00 

cm 

0 

.00 

cm 

(0.00 

in) 

0 

.00 

cm 

(0.00 

in) 

0 

.00 

i n 

(0.00 

in) 


o c 


o c 


1.27 

cm 

(0.50 

in) 

38. 10 

cm 

(15.00 

in) 

12.70 

cm 

(5.00 

in) 


Swept Wing: 
sweep 

quarter chord MAC location: 
fuselage station 
water line 
butt 1 ine 

Straight and Swept Wings: 

wing: 

area 

aspect ratio (on nominal chord) 
span 

nominal chord (constant) 
quarter chord water line 
twist 


25° 

6.6*4 cm 
38.10 cm 
12.70 cm 


517.00 cm- 
5.00 

50.80 cm 
10.16 cm 
38.10 cm 
0 ° 


25° 

(2.71 in) 

(15.00 in) 
(5.00 in) 


(0.556 ft : ) 

(20.00 in) 

( *1 . 00 in) 

(15.00 in) 
0 ° 


ter line 0.0 is snail tunnel floor w nth node l on tunnel 


renter line. 


3P 


Con t i nued 






TABLE I ~ 

MODEL DIMENSIONS 


wing and tips: 
area 

aspect ratio (on mominal chord) 
span 

nominal chord 
leading edge slat: 

area (projected onto maximum chord): 
w i ng on 1 y 
w 1 ng and tips 
span : 

wing only 
wing and tips 
chord (maximum) 
slot width 
def lection 

trailing edge flap: _ . 

area (projected onto maximum chord) 

span 

chord (maximum) 
slot width 

deflections (wing chord to flap 
upper surface) 
straight wing 
swept wing 


968.00 cm : 

6.00 

76.20 cm 
12.70 cm 


103.00 cm 2 

155.00 cm 2 

50.80 cm 
76.20 cm 
2.03 cm 
0.127 cm 
80.00 

23*4-00 cm 2 
50.80 cm 
*4 . 60 cm 
0.0*4) cm 


76 . 00 ° 

60 . 00 ° 


(1.0*42 ft 2 ) 

(30.00 in) 
(5-00 in) 


(0.1'i Ft 2 ) 
(0.16/ ft 2 ) 

(20.00 in) 
(30.00 in) 
(0.80 in) 
(0.050 ir>) 


(0.252 ft 2 ) 
(20.00 in) 
(1.8) in) 
(0.016 in) 


76.00° 

60 . 00 ° 




Figure 2.1 influence matrices for sou 















igure 2.2 Application of influence matrices 

























FITTED (ITERATIVE METHOD) 
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Performance of 
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in 164 • X 23fc ' tunnel). 
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(a) AXIAL VELOCITY DUE TO VORTICES 



x STRAIGHT WING O * ?S" SWEEP, > - 0° 
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Figure 2.6 Effect of combined 'weep and anqle-of -attacL 
on roof-mi nus -f loor influence curves 
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Figure 2.7 Sensitivity of MATCH" and iterative methods 

to data scatter ("SA$ , ‘ car model in 164* x 234' tunnel) 
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Figure 2.9 Application of the "Least-Squares" Approach to "noisy" Data 
(SAS Car Model in l6i' x 23£’ Tunnel) 
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SLOT AND PLENUM: 
TYPICAL CROSS SECTION 




Figure 3.7 Wal. ..ressure orifice locations for experiments with 
the knee-blown flap models. 
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SWEPT, KNEE-BLOWN 
JET-FLAP MODEL 
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Figure 5.2 Correction of lift and drag data for the swept, 
knee-blown, jet-flap model. 

a) Sensitivity to sweep of line singularities 
("cross" effects included) 
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Figure 5*2 (continued) Correction of lift and drag data for 
swept * knee-blown, jet-flap model 
b) Sensitivity to 'cross' effects. 

(Swept singularities used set at 15° angle of attack) 
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Upwash at P, due to A, B, C etc, is neutralized by 
downwash due to A', B', C* etc. 



Sidewash at P, due to R, S etc, is balanced by 
sidewash due to R' , S' etc. 


Figure 6.1 Occurrence of 'neutral 5 points in 
a wind tunnel cross section. 
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APPENDIX - I 


i 
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PROGRAM DESCRIPTION 

Capabilities : 

The tunne 1 ~wa 11 -ef f ec t correction program is a generalized version to 
handle complex pressure signatures arising from powered model tests. It 
essentially solves an inverse problem of determining the strengths of 
potential singularities, the geometry of which has been specified, to satisfy 
the measured pressure signatures on the tunnel boundaries. The number of 
s i ngu I a r i t i es can be fewer than the number of pressure signature points since 
the present approach satisfies the boundary condition in least squares sense. 

It is possible for the user to specify arbitrary orientations and 
geometry for the potential singularities to model the actual flow as closely 
as possible. In the present version of the program, no assumptions regarding 
the symmetry or an t i -symme t ry of the influence coefficients are made to re- 
solve the signature into vor tex- re I a ted and source-s i nk-re la ted parts. This 
resolution is done iteratively during the numerical computations. At 
present, the tunnel geometry however is restricted to rectangular shapes, 
since the computational procedure uses imaging technique to ensure zero- 
nj»mal-flow through the tunnel walls. However, alternative arrangements are 
made for cases involving the *4(3 1 x 8 0 ' tunnel cross section (see below). 

Normally, the difference between the observed supervelocity on the roof 
and the one on the floor is used as boundary condition for obtaining the 
vortex strengths. However, with powered models involving jet impingement 
on the floor, it may not be desirable to use the floor signature in the 
calculations. In such cases a flag, IRE, car be set to handle only the 
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roof signature. Note that this requires that the cross term flags (KROSG 
and KROSQ) should be turned on and that the number of iterations (ITERMAX) 
should be set larger than unity even if all singularities are placed 
symmetrically with respect to the tunnel cross-section. 

The program coding was developed using a VAX- 1 6 computer. FORTRAN 
statements that may cause problems in other systems are identified by the 
characters VAX in columns 73 ~ 75 of those statements. When using other 
systems these statements should be appropr i a tel y replaced. 

The present coding is written with the assumption that the pressure 
signatures and load coefficients to be corrected are made available in a mass 
storage file. The subroutine READER reads in these values using FORT RAM I/O 
unit number 10. This subroutine is written to handle specifically the 
KBF model data of Lockheed-Georg ia . In this case, eight rails of tunnel 
wall signature data were available in a mass storage file in the form of 
super-veloc i t ies rather than Cp-values. Also, since the x-wise locations 
of pressure points for rail No. 7 was different from the rest of the rails, 
a subroutine INTER is employed to linearly interpolate the rail-? data to 
the standard x-wise locations. Since the general user’s data structure 
will be different from that of Lockheed's KBF tests, these subroutines 
and their calling sequences in the driver program might have to be replaced. 

Preparation of Input Data 

The overall sequence of computations and the effects of different flag 
settings are shown in the flow chart given In figure A1 . The meaning of 
all input variables are explained in the next section. A typical run of the 
program involves one of the two cases: (1) The required matrices are all 
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Figure A1 . 


Flow chart for wa I 1 -pressu re signature-based 
tunnel interference program, 
a) Pre-analysis routines 
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available and only the signatures need be processed. (2) The matrices must 
be generated and saved for future use. The input sequence for these two 
cases are as Follows. 0) Matrices available: Prepare Card Number 1 
through 3 as indicated under next section. Set MATSAV-2 in Card Number 3- 
Skip Cards 4 through 9 and prepare Card No. 10. (2) Matrices must be 
generated: Prepare all cards. No. 1 through No. 10. Set MATSAV-' in Card 

No. 3- 

Description of Input Variables 

Main InpiU - The main input portion consists of a title, all flag 
variables and a few key variables related to the model geometry as outlined 


below. 


1 I T! 


Format ““ 80Al 


TITLE: 


Test Description 


2 ITERMAX MATPRT MATSAV I PRT KROSG KR05Q, 1C0RR JETEFCT .ormat 


ITERMAX: 


Maximum 


number of iterations to be performed when cross- 


effect- terms are to be included in analyzing the signa 
tures. Program automatically sets this to unity if both 
KROSG and KROSQ are zeros. 


MATPRT: 


A non-zero value causes all influence coefficient matrices 
to be printed out. 


MATSAV : 


Three-wjy flag. 

-0: Generates matrices but does not save them. 

-1: Generates matrices and writes them to FORTRAN unit 

Number 8. 
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IPRT: 


KROSG : 


*2: Implies that all required matrices are available 

from a previous run. They will be read from 
FORTRAN unit No . 8 . 

Flag for printout detail. 

“0: Prints out a one page summary for each tunnel test. 

“1: Prints out details for each iteration. 

When non-zero, calculates and takes out the sidewall 
upwash due to vortices in determining the u-velocity 
boundary conditions on side walls from pressure 
coef f icients . 


KROSQ: When non-zero takes out the cross effect of sources/ 

sinks from the Roof/Floor signature* before calculating 
the circulation strengths of the vert ires, (Note: 

A non-zero value is meaningful only if sources/sinks 
are not placed synVDe t r i ca 1 I y with respect to the tunnel 
cross sect ion) . 


I CORR : The centerline Inter, arence velocities at the x-location 

correspond i ng to this index will be used In making 
final corrections to angle of attack and the loads. 

JETEFCT: Non-zero if model includes a lifting jet. 

(See Part II of the report) 


SAREA AWB BWB Forma t — 8F10 . ^ 

SAREA: Model Reference Area used In normalizing the load 

coef f Icients. 
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Constants in the But l er~Wl 1 II ams Equation for correction 

BWB : , , v 
to angle of attack, Aa » AWB-CL c /(l + BWB ;; Cp C ) . 

Geometry Input: The input of this section pertains to the tunnel geometry 

and the singularity geometry. 


This entire section should be skipped while preparing the input in MATSA\/=2 in 
the Main Input which implies that all requ i red matrices are already available. 


LAYERS 

IRF 

NR 

NW 

NV 

NS 


LAVERS: 

IRF: 


NR: 

NW : 

NV: 

NS: 


XPVOR: 


No. of Image Layers to be used ( Re commend ed : 5) 

Flag for determining whether floor signature is to be 
used or not. 

=0: Implies usage of roof signature only. 

= 1 : Implies usage of (roof-floor) values. 

(Note: If roof signature only is to be used, set KR0SQ= 

and ITERMAX-1) 

No. of roof signature points. 

(it is assumed No. of floor signature points is same) 

No. of side-wall signature points. (Both sides are 
assumed to have some no. of points). 

No. of vortex singularities 
No. of source/sink singularities 


5 

B 

1 H 

XPVOR 

XPSRC 


B: 



Tunne 1 


H: 



Tunne 1 


Pivot point for pitching swept vortex, normalised w.r.t 
Span (See sketch). Meaningful only when the vortex is 
swept . 
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Pivot point for pitching swept source/sink normalized 


with respect to semi-span, definition is similar to 
XPVOR given above* 



F o rma t - - 8 F i 0 * ^ 



Forma t --8F 10.^ 


2 R 


NR 


ZL, 


21 . 


Format-'-8F10. ^ 

Forma t --8F 10.^ 
Forma t --8F 1 0 . h 


2L W ) Format--8FI0. *4 


Non-dimensional coordinates, 
x/8 , y/B, 2/N, of the roof / 
floor signature points. 


Non-d imens I ona i coord i nates , 
x/B f y/B, z/K of the side-wall 
signature points. 












8.1 

8.2 

8.NV 


XV, 

I 

YV 1 

ZV , 

SBV j 

PSIV, 

ALFV, 

XV 2 ! 

YV 2 

ZV 2 

sbv 2 

PS 1 v 2 

* LF v 2 j 


Forma t~-8F 1 0 . 4 


XV,. 

YV.„. 

ZV, (W 

SBV. 11( 

PSIV, I4 , 

ALFV..,. 

NV 

1 nv ! 

NV 

NV 

NV 

NV 


Forma t~-8F 10-4 


(XV,YV,ZV): Non-dimensional coordinate x/B, y/B , z/H of the "roof” 

point of horse-shoe vortices. 

SBV: Vortex span normalized with B 

PS IV: Sweep angle for the vortex (degrees) 

ALFV: Pitch angle for the vortex (degrees) 


9.1 

9.2 
9. NS 


XS . 

YS , 

ZS 1 

5BS, 

PSIS j 

ALFS ] 

XS 2 

YS 2 

ZS 2 

sbs 2 

PSIS 2 

ALF5 2 


Forma t — 8 F 10.4 
Forma t--8F 10.'+ 


XS NS 

ys ns 

zs ns 

sbs ns 

PS,S NS 

alfs ns 


Forma t--8F 10.4 


The definitions of these variables, defining the source/ 
sink locations and geometries, are similar to the ones 
for vortices. 

Test Signature input : This last input card contains the values of key 

variables identifying the test, the signatures correspond ing to which are 
to be picked up from a mass storage file. This card can be repeated as many 
times as desired to process all required runs. (NOTE: The u$cr need to 

replace this section of the coding. See earl ler comments about subroutines 
READCP and INTER) 
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10 


(TEST 

IRUN 

1 PH IN j 

i PMAX 

1 FLOOR 

1 ROOF 

IWALl 

IWAL2 


1615 


I TEST : 
1 RUM : 

I PM IN: 
1PMAX: 


! FLOOR 
I ROOF 
iWALl 
IWAL2 


Test Number 
Run Number 
Point No. , minimum 
Point No.* maximum 

For the given run number the program will process data for 
all points i in the range 1PM1N < i < ? PMAX . 

Pressure signature rail no. for floor 
Pressure signature rail no. for roof 
Pressure signature rail no. for first sidewall. 

Pressure signature rail no. for second sidewall. 


Mass-Storage File Requirements: 

In addition to the standard input/output FORTRAN units (#5 and if 6 in the 
coding)* the coding employs four other mass-storage files* as explained 
below. 

UNJT-7: Output file. The test no., run no., point no. are written 

to this file along with a summary of measured and 
corrected angle of attacks and load coefficients. It 
may be used in preparing plots if so desired. 

UNIT-8: Input/Output file. 

This would contain all the input data entered in the 
section “geometr/ Input 11 , all the required influence 
coefft. matrices and the least square i nverse mat r i ces . 
This file has to be generated and saved when the program 
Is run for the first time or whenever a change In any of 
the input variable described in the section “Geometry 
lnput M is made. 


106 



UNIT-9: 


Output file 


For each data point and for each iteration this would 
contain the input signatures, the recalculated boundary 
conditions and the centerline interference velocity 
components. It may be used in preparing machine plots 
to evaluate the program. 

UN I T- 1 0 : I nput file. 

For each data point, this file should have the pressure 
signature. The structure of this file is left to the 
user (See comments about subroutine READCP above). 

Output Format 

A complete sample output follows the program listing. The output is 
sufficiently well annotated for easy comprehension of the print out. In the 
printer plots of the input signature, calculated wall supervelocities and 
tunnel center line velocities, the cor respondence is readily established by 
looking for the same plot symbol under the tabulated data. In the table of 
corrected a and load coefficients, the values labelled "CLASSICAL" are the 
ones obtained using Butl er-Wi 11 i ams equation. In the output annotations the 
word ROOF would mean either roof alone or (roof-floor) values depending upon 
how the input was arranged. Notations like "U-Q-CL" imply "u-velocity due 


to sources at tunnel center line." 
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EXAMPLE. . XBF_TESTS . . .USE AVAILABLE MATR ICES .. EXCLUDE CROSS TERMS 
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3 CONTINUE ‘ CRDTEG 045 

RETURN CRDTEG 046 

END CRDTEG 047 
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CONFIGURATION: 
FLAP BLOWING: 
ANGLES-OF-ATTACK.: 


i 

-0.4 


r» . 

OF POOH ; 

APPENDIX 3 

Examples of Signature Analysis 

STRAIGHT-WINGED KBF MODEL, NO TIPS. 

C v = 2.0 

6°, 12°, 18°, 2A° 

PLOTTED AFTER THIRD (l.E. FINAL) ITERATION 




I 


• • © $ 

® 0 

• 



0.0 C . 4 

X/B 

C . L 


■ • - 

1.6 


151 


a = 6 


TEST 62, RUN 45, POINT b 

U,SI0E Wftll: <*> MEASURED, © C3LC JLfi'EQ U. CENTER 

0.1600 — r - - - 


0.1200 


o 



U, ROOF-FLOOR: <j> MEASURED. © COLC'JLO’ED U, CENTER 

0.6000 , - - - ' ' 



- 0.2000 




- 0.0400 

U. ROOF-FLOOR: <j> MEASURED , 0 CALCULATED W, CENTER 

0.6000 - - 



- 0.2000 
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APPENDIX 4 


LEAST SQUARES APPROACH FOR THE NASA 40 ' X 80* TUNNEL 

The image method employed by the LSQITER program cannot be applied 
directly to the NASA 80 1 x 40 ‘ tunnel because of its ncn-rectangu! ar cross 
section. The influence coefficients required in the LSQITER program have 
to be generated using an alternate approach. This Appendix presents the 
results for interference factors for the 40 ‘ x 80 1 tunnel and explains how 
these are used to construct the influence coefficient matrices required by 
the LSQITER program. 

Influence factors due to an isolated singular! ty 

The influence factors due to a single horse shoe vortex or a finite 
length line source are obtained by using the vortex panel method described 
in Ref, 9. Figure A4.1 shows the theoretical flow model used to generate 
these factors. A length of 288 * of the funnel is panelled with vortex 
latices. Velocities due to a centrally located singularity are calculated 
at these panels. Panel circulation strengths are then obtained which 
satisfy the zero normal velocity condition at the tunnel surface. 

Tunnel wall super velocities are then computed as the sum of panel- 
induced and s i ngul a r i ty- i nduced effects. These calculations are done at 
various values of x/8 at the roof, floor and the sidewall locations indicated 
in figure A4 . 1 . The center-line influence factors are computed by omitting 
the effects of the central singularity and including only the panel circu- 
lation effects. The supe rve loc i t i es thus computed are normalized by the 
factor Q/C for cases involving sources and by the factor 2Pb/C for cases 
involving horse-shoe vortices, where C is the tunnel cross section area and 
b is the singularity span. The results for sources and horse shoe vortices 
of different *pans and for both horizontal and vertical orientations are 
presented in Tables A4 . 1 through A4 . 6 . 
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Generation of influence coefficient matrices 


Using the normalized influence factors presented in the tables, influ- 
ence coefficient matrices are constructed. The elements of these matrices 
are of the form a; j where a is the induced velocity due to j-th singularity 
(of unit strength) at i-th point. Thus the influence factors given in the 

table must be multiplied by l/C for sources and by 2b/C for vortices for 
ongoing use. 

The independent variable, x/B. presented in the tables corresponds to 
stream-w.se locations on a local coordinate system whose origin is at the 
singular. ty. Once the pressure ports locations are chosen and the positions 
of singularities at the tunnel centerline have been selected, the relative 
streamwise distance between a given singularity ( j ) and the pressure measure- 
ment point (.) is known. This relative distance normalized by the tunnel 
breadth, B, is used as the independent variable to pick values from the 
tables. It may be necessary to interpolate the tabulated values. 

The LSQITER program has been written to handle the most general cases 
involving singularities that could be swept, pitched and be located off- 
center in the tunnel. However, due to restrictions on time and effort the 
influence factors for the kO ' x 80 ' tunnel have been generated only f or 
cases where the singularities are unswept and are placed midway between 
roof and floor. Consequently, many of the influence coefficient matrices 
required by the LSQITER program become null matrices. The following list 
defines the matrices required. 


UGRF : 

u 

due 

to 

r. 

UGWL : 

u 

due 

to 

r 

* » 

WGWL : 

w 

due 

to 

r. 

UQWL: 

u 

due 

to 

Q, 

UQRF: 

u 

due 

to 

Q. 

UQCL: 

u 

due 

to 

Q, 

WQCL : 

w 

due 

to 

Q, 

UGCL: 

u 

due 

to 

r. 

WGCL: 

w 

due 

to 

r. 


Non-zero 
Ze ro 
Ze ro 

Non-ze ro 
Zero 

Non-zero 

Zero 

2ero 

Non-^ero 
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Note, however, that the nul! matrices must be made available to the LSQITER 
program with all elements set to zero 

Data file structure for the influence coefficient matrices 

The influence coefficient matrices generated for special cases like the 
40’ x 80 ' tunnel must be made available to the LSQITER program via FORTRAN 
input UNIT NO. 16. The structure of this data file is as follows: 

The first few lines of data correspond to what was described as 
Geometry Input' in the input description of the program given in detail on 
pages 103 through 105- Input line numbers 4 through 9 must be defined in 
this data file accordingly. The FORMAT s for the variables are 1615 for 
integers and 5E16.8 for real numbers (see also subroutine TAPEI0 in program 
listing, page 132). Note that the variable "LAYERS" loses its significance 
and that there can be no sweep or pitching of the singularities. Following 
these lines of input, the data file must now contain the elemenets of the 
matrices listed above in the same order. The elements (a:j) of each matrix 
must be sequenced such that the subscript i varies more rapidly than the 
subscript j. (See program listing on page 132). 

Running the LSQITER program for the 40' x 30' tunnel 

Once the data file containing the influence coefficient matrices is 
constructed as described above, the LSQITER program can be run to process 
the 40 1 x 80' tunnel signatures. The input on UNIT 5 is identical to the 
rectangular tunnel case input described in pp. 101-106 with the following 
except i ons : 

1. The program must be signalled to expect the special influence coefficient 
matrices. This is done by assigning a value of 3 to the flag MATSAV in 
input card number 2 described on page 101. This causes the program to 
read these matrices from FORTRAN UNIT 16 instead of calculating them 
through imaging techniques. 

2. At present the flags KR0SG and KR0SQ in Card-2 must be set to zero since 
influence coefficients for cross effect terms are not available. 
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«.„ the Geometry -put Section (Card-b through Card-S) at these will 
now be teed from the matrix data file via UNIT lb. 

addition to the mats stora 9 e files described on paries 106 and 107, 
the matrix data file described in this appendix must be pre-ass,gne to 
FORTRAN UNIT , 6 . For a given geometry of singularities and wa 1 pre 
points, the least sguare inversion process needs to be done only once 
program writes on, ... matrices on UN1T-8 which has to be saved for fn ure 
OS, Subsequent runs can be made with NATSAV-2 and the spec, -flue 
coefficient matrix data file need no. be made available via UN.T-,6 (See 
comments on mass stroage file. UNIT-8 on p„g- '081. 
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TABLE AV!B Model-center Interference coefficients for horizontal line 
sources in the kO' k 80 1 tunnel. 
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TABLE A4/2A Tunnel roof influence coefficients for horizontal horseshoe 
vortices in the 40 1 x 80 1 tunnel. 
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Table shows interference Aw/U , at model centerline, for T/OJ^C/b) - 0.50 
(Values are the same as conven t i ona 1 1 y -def i ned 6) 
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TABLE AA/3A Tunnel roof influence coefficients for vertical line sources 
in the bO 1 x 80 1 tunnel. 


SOURCE SPAMS 
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’ABLE AA/3B Model-center Interference coefficients for vertical line 
sources in the ^0 1 x 80 ‘ tunnel. 


BOUND VORTEX SPANS 
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TABLE A4AA Tunnel sidewall influence coefficients for vertical horseshoe 
vortices in the b0‘ x 80 1 tunnel. 
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TABLE A4/^B Model-center interference coefficients for vertical horseshoe 
vortices in the *+0* x 80 1 tunnel. 
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